Complex base numeral systems. 



Jarek Duda 



Jagiellonian University, Reymonta 4-, 30-059 Krakow, Poland, 
email: dudaj@interia.pl 



Abstract 



In this paper will be introduced large, probably complete family of com- 
plex base systems, which are 'proper' - for each point of the space there is a 
representation which is unique for all but some zero measure set. The con- 
dition defining this family is the periodicity - we get periodic covering of the 
plane by fractals in hexagonal-type structure, what can be used for example 
in image compression. There will be introduced full methodology of analyz- 
ing and using this approach - both for the integer part: periodic lattice and 
the fractional: attractor of some IFS, for which the convex hull or properties 
like dimension of the boundary can be found analytically. There will be also 
shown how to generalize this approach to higher dimensions and found some 
proper systems in dimension 3. 

1 Introduction 

Standard numeral systems allows us to represent points from M in some base (z) as 
a sequence of digits: 



In this paper will be shown how to generalize it to 2 dimensions, using the structure 
of complex plane C = M + jM. It was previously done by Donald Knuth [Tj , who 
has introduced imaginary base systems, or in |i2j| where the j — 1 base system was 
introduced and where was shown that such representations can be used to simplify 
operations involving complex numbers in today's microprocessors. 
In this paper will be presented much larger family of complex base system, which 
seems to be complete. 

We expect from the numeral systems, that the representation function is surjec- 
tive - each point has at least one representation. 

But we rather cannot expect injectivity - for example 0.111111. ..2 = 1.000000. ..2 - 
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there is countable number of points that have two different representations. 

In the 2- dimensional case such set won't be countable - it will have Hausdorf's 

dimension in [1, 2) range. I will assume only, that it's Lebesgue's measure is 0. 

Definition 1. We will call the numeral system proper, if the representation function 

is surjective and pseudoijective: 

^Scc '■ IJ'2{S) = 0,\/xec\s X has at most one representation. 

Choose some N 9 n > 2. 
We will use n := {0, .., n — 1} digits and some base z € C. 
We have to analyze two completely different sets: 
Integral part: I^ = {Yli>o^i^^ : ai e n, 3ArVi>jv = 0} 
Fractional part: — {X]j<o ^^-^^ : G n} 

For the binary system n — 2, z — 2, fi — {0, 1}, I2 — I', F2 — [0, 1]. 




Figure 1: Integer part we can get using 6 youngest digits for n = 2 and z — j + 1 
(left) or z = j — 1 (right). It looks like in the first case we will make 'spirals' around 
0, but in the second we should cover the whole space . 

We will check, that they are fulfilling self similarity equations: 

z-F,^ [j {F, + i) (1) 

i=0..n-l 

U {zh + i) = h. (2) 

«=0..n-l 

For the binary system 2[0, 1] = [0, 1] U ([0, 1] + 1), 2Z U (2Z + 1) = Z. 
We will check that the necessity conditions implies: 
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n=2, D=-2 



n=2, D=-1 



n=3, D=-2 



Figure 2: Some fractional parts. The color represents two oldest digits. 
Lemma 2. If [n, z) is proper 

/X2(F.) > 0, /i2(F, n (F, + 1)) = 0, 

n = \z\'^. 

It's why we will be dropping the upper index of F". 



(3) 



So Fz is some compact, positive measure, central symmetric ([3]) set and when 
we take copies translated by all points of some discrete set Iz we 'tile' the whole 
plane, such that tiles should intersects with its neighbors only on the boundaries. 
Intuition says that it should be periodic tiling: 



(4) 



We will see that it means, that z = D/2 + j^/n — [D/2Y 
for some D eZr\ 

Hypothesis 3. All proper (2- dimensional) numeral systems are periodic Q). 

I couldn't prove it. On Fig. [3] are numerically found Hausdorff's dimensions of 
F^^j^ for Lp G [7r/2,7r]. The tree peaks with dimension two are the periodic cases, 
while a — >^ TT, -F^ is becoming a segment - its dimension drops to 1. 



We will concentrate on the periodic case in this paper, and show that the picture 
looks like it seems (like on Fig. |2] or Escher's pictures) - we have slanting two 
dimensional lattice {X = Z + zZ,) and each of its parallelogram (build of 1 and z) 
corresponds to exactly one copy of F^ (so [Ji2{.Fz) = \z\ sm{ip)), which is connected, 
simply connected, F^ = int{Fz) and it intersects with its 6 neighbors (for if ^ ±7r/2) 
on its boundary, which dimension can be calculated analytically. 
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Figure 3: Numerically found Hausdorff's dimension of F^. 



These boundaries (SFz) are the sets on which points has two or three representations 
- we can define dense set on which we are loosing injectivity - like for the binary 
system: 

S=\Jz'' ih + SF,) 

In one dimensional case, this set is countable - has dimension equal 0, it's less by 1 
than space dimension. The skewness makes that this time, we have a bit worse (see 
table on page 16). 

The surjectiveness condition can be written in a few ways: 

C = F, + I, = J2^'F, = Y^ (5) 

i>0 i<0 

we can check it: 

• by checking that Iz — 'L^zL 

• by checking that e int(F^) 

In the next section will be introduced formalism and proven some basic 
properties: selfsimilarity relations, \z\^ — n condition, that it's enough to analyze 
argument of z in one quarter and be characterizedS the periodic cases. 
In the third section we will focus on the integer part. We will show in what cases 
= Z + 2;Z, what is sufficient for system to be proper and that in the rest of cases 
situation isn't much worse. The main tool will be the reduction - the analogue of 
shift right with throwing out the youngest digit. 

In the forth section we will focus on the fractional part. Firstly there will be 
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presented results from [3] - a methodology of finding the convex hull of and 
its properties analytically. In the second - main subsection, will be introduced 
the methodology of constructing the approximation of the boundary of F^, which 
allows to show that the tiling looks 'nice', quickly approximate F^ or to calculate 
analytically the dimension of the boundary. 

In the last section, will be briefly presented a proposition of generalization of this ap- 
proaches into higher dimensions and shown a family of proper system in dimension 3. 

We can use the complex base systems is image compression too - in hexagonal- 
type structure, we should get better correlations. 

The other thing is that calculating the transform of the whole image would cost 
too much, so the image is usually split into a periodic lattice of squares, which are 
processed separately. In this approach lossy compression creates recognizable lattice 
of lines. Splitting into lattice of fractal-like shapes shouldn't have this problem. 
To encode inside such shape - {J2o<i<N ^i^^ '■ ^ '^^^ wavelet transform 

as it was just usual binary system, for example basic Haar wavelet for n = 2 would 
be: fk = 2afc - 1 if Vi>fc = 0, else. 

The boundary of larger integer sets is created from boundaries of smaller ones, so 
the wavelet functions vanishing on the boundaries should behave well in this case. 



2 Basic definitions and properties 

Definition 4. Fix some 2 < n eN, z E C : \z\ > 1 

Digits n := {0, 1, .., n — 1}, 

Pseudorepresentatnions P := {(ai)i6z : VjOj E Z, 37vVj>iv CLi = 0}, 

Shift left(right): shr(aj) := (ctj-i), shl(aj) := (aj+i), 

Representations P" := {(cti) E P : Vjaj E n}, 

P+ := {(a,) E P : V,<o = 0}, P"+ := P+ n P", 

P- := {(a,) E P : V,>o a, = 0}, P"" := P" H P", 

Representing function P3 (oj) {cii)z '■= 

Integer part (/^ =)-^" '■= Pz^, 

Fractional part {F^ =)P" = P"^- 

We will use standard set arithmetic: 
A + B = {a + b:aEA, bEB}, f{A) = {/(a) : a E A}. 

Observation 5. .■ 



(shl(aj))^ = z ■ (ai) 



(shr(ai))^ 



= z 



-1 



shl {P"-) = P"~ + (..,0, n, 0, ..), shl (P"+) + (.., 0, n, 0, ..) = P 
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So we have the selfsimilary equations: 

z-f: = f^+% z-i:+n = i:. 



We've just checked Q and ^ conditions. 

The surjectiveness condition means that P" = C 

P^ = Y^z^F, = F, + h (6) 

A;>0 

it's the sum of countable many translated copies of F^ so we must have fi2{Fz) > 0. 
Pseudoinjectivity implies that Vj=i..„_i/i2(Pz H (P^ + «)) = - points in this intersec- 
tion would have at least two representations - with and with i at zeroth position. 
Multiplication by z is rotation and rescaling. Rotation and translation doesn't 
change the Lebegue's measure, and the scaling multiply it by the scaling factor to 
power dimension - from the self similarity condition we have: 



U (p^+A 

\i=0..n-l ) 



n ■ fi2{F,) = /i2 I [ J (F, + 2) I = fi2{z ■ F,) = |z|>2(P. 
So we've checked Lemma lU 

Observe that if we assume \z\'^ = n, from the derivation above: iX2{Fz) > auto- 
matically implies pseudoinjectivity - we have only to check the surjectiveness. 
For the rest of the paper denote: 

z = V^e^"^. (7) 

When we conjugate z, F^ will be conjugated. 
When we take —z, we get {—n = n — [n — 1)) 

= E.<o = E.<o - E... .<o An - 1) = - = Fz- 



(n— 1)2 ^ p {n—l)z 



Observation 6. 

F, = Fz (8) 

^ ^ {n — l)z , , 

F.z = Fz- (9) 

We will now find possible z, fulfilling periodicity condition {z"^ & 7j + zZ) and 
\z\ = n. 

Thanks of ([s]), we can restrict to positive imaginary part and eventually take con- 
jugation of everything. 



d + i \Jn — (P' for some d G 



2; = d + J 



„2 



z 



d"^ -n + d"^ + 2jd\/n - d^ 



Comparing imaginary parts, we have D := 2d E we finally have: 
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Observation 7. All periodic cases are: 

z = D/2+jy/n-{D/2y (10) 

= Dz-n (11) 
for some 2<nGN, D eZf] {-2^, 2^) . 

In the next section we will prove: 

Theorem 8. Complex system fulfilling the periodicity condition is proper iff D < 1. 

3 The integral part 

In this part we will check whether = X, where we will denote 

X ■.= 'L + z'L. 

Because Iz + Fz = IJfc<o "^^^z-, we'll have checked the Theorem [s] then. 

3.1 The general behavior 

Fix some n, z as in Observation [7j 

We have n > 2, D"^ < An - it's easy to check: 

D<n, n = D ^ n G {2, 3} (12) 

The equation z"^ = Dz — n can be written as 

(n,-Z},l). = (13) 

where we use the natural identification: 

Zf'^ 3 {bo, ..hk-i) = (dj) : Vj<ovi>fc = 0, Vo<i<fc Oj = bi. 

Now take some C 9 x := (6o, ■■■,bk)z, where Vj6j G Z. 

If we add to this pseudorepresentation of x any multiplicity of (n, —D, 1) or some of 
its shift, it will still represent x. If we subtract [bo/nl {n, —D, 1), we will make the 
zeroth digit to be in n. 

It suggest an algorithm to change a psudorepresentation into representation - use 
above construction on succeeding digits to bring them into n. 
The question is: will it finally stop? 
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Example 9. We can use (n, —D, 1)^ = in many ways. 

n = 3,D = 3^z=^+jy^, (3,-3,l), = 

We can change some pseudorepresentation into representation: 

(7, -8, 7, -2) (1, -2, 5, -2) ^ (1, 1, 2, -1) ^ ... ^ (1, 1, 0, 0, 2, -1) 
Find the representation of some x G X: 

x = 10z-12 = (-12, 10)^ = (0, -2, 4)^ = (0, 1, 1, l), = z + z^ + z^ 
Or find the position on X for some representation: 

z^ = (0, 0, 0, 0, 1), = (0, 0, -3, 3), = (0, -9, 6), = (-18, 9), = 9z- 18, 



Digression We can use arithmetic in such bases similar as usual - for example 
adding would be adding two representations position by position - we get a 
pseudorepresentation of the sum. We have to change it into representation. In the 
binary system we would use (2, —1)2 = 0, this time (n, —D, 1)^ = 0. 
We can create parallel algorithm for it - in the binary system we divide the 
representation into two parts and sum them parallelly, for every possible carried 
digit (0 or 1). The difference in this case, is that we have more 'carried cases': for 
addition from 8 (for Z) = 1) to 21 (for n = 3, D = ?>). 

X E 2, + zIj = Xq + zxi can be 'psuedorep resented' as x = (xq, Xi)^, so 
Lemma 10. For any x G X, there exists its reduction r(x) := y E X , such that: 

X — yz E n. 

We can think about the reduction as throwing out the youngest digit (g n). 

We will now check if r have a nonzero fixed point (s). 
Assume that r(x) = x for some x = a + zb: 

r ((a, b, 0)z) = r ((a - kn, b + kD, -A;)^) = (6 + kD, -k)^ 
where k = \_a/n\. 

a + bz = {b + kD) — kz, so b = —k, a and b has different signs, 
a = b + kD = b — bD = 6(1 — D) - we must have D > 1, 



b = -lb{l- D)/n\ = \b^], but < ^ < 1 (12) so 6 > 
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for 6 = 1 we have (1 — D, 1) for any D >2, 

for 6 = 2, we would need 2{D — 1) > n - it can happen only for n = 3, D = 3 - it's 
very special case we would focus later on. 



Because of (12) it's the largest possible b. 



We have that there are 1-3 fixed points of r. We would like to show that from 
all points of X we would get in finite number of reductions to one of these points. 



Lemma 11. |x| > ^/n + 1 =^ \i"{x)\ < \x\ 

Proof- \r(x)\ = l^^l < — kl+(V"+i)(v^-i) 



< 



\X\. 



□ 



Because we are in the discrete lattice X, starting from any point x G X after a 
finite number of reductions, we will get into the ball 5 := {x G X : |x| < ^Jn + 1}. 
It's finite set, all fixed points are in it - if we would show that there are no cycles in 
it, the reduction process will always stabilize. 

Definition 12. Attractor of fixed point s, is A^ := {x : 3jr*(a:) = s}. 

B is finite, so for a finite number of cases we can do it by checking that in fact 
all points of B will be in some attractor. I've done it for n = 2, 3 (Fig. 111). 
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Figure 4: All cases for n = 2,3. Different colors represents different attractors. 
Black dot is the center of coordinates. 



For n > 4, ^/n + 1 < n — 1, so?7, + zn, which is in the attractor of 0, covers 
{x E B : arg(a;) G [O,^?]} (see Fig. 5). 
V? G (0, vr) : if arg(a;) G [n,2n), arg(r(a;)) < arg(a;). 
So finally this cycle could be only in {x G -B : arg(x) G {ip, vr)}. 
li (f > n/2, {D < 0), after one reduction we get into the n + zn, else after 2 
reductions we get into the second stationary point. 



We've just proven Theorem [81 



□ 
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3.2 Special cases 

We will now concentrate on the D >2 cases: when ^ X. 

In this case we can split X into two or three disjoined subsets (attractors of fixed 
points) - we can still use is as a numeral system, but this time we have to remember 
the number of fixed point too (or of cycle for some cases in higher dimensions). 
For each fixed point s = r(s), we can split the attractor into subsets (Ak) of points 
which have the same number of reduction [k) to get to the fixed point (we need k 
digits to encode it): 

zA'; + n = A':+\ r{A':+') = r{A,) 

So because A^ = {s}, 

:= sz'' + z^n (= sz'' + A'^) 

0<i<k 

Using As = J2k>o "^s C A] C ... C C we have 

zAs + n = A,, (14) 
it's the same similarity equation as for ([2]). 

There can be some nice symmetries here: check for which xq, xq — Iz fulfills (§ 
z{xo — Iz) + n = z{xq — Iz)+n— 1 — n = zxq + n — 1 — {zlz + n) = zxq + n — 1 — Iz 
If we want the right side to be equal xo — Iz, we must have zxq + n — 1 = xq 

xo = J — - (15) 
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Figure 6: The difference for D = ±3. 



For D — 2,we have xq = , " = j\/n — 1 — z — 1 e X, 

For n = 3, = 3, = 22; — 4 we can say, that A^-^ separates this two copies. 

Unfortunately for the rest cases ^ ^ - we can see it using 

Al= '^^= An-l-n)= z\n-l)-Al = {n-lf^-Al 

i=Q..k-l j=0..fc-l i=0..fc-l 

so A^g = sz^ + (n — 1)^17 — Aq, but this point moves with k in the other cases. 
Finally 

Theorem 13. We can have the following periodic cases: 

• For D < 1, X = Aq, 

• ior D ^2, X ^ AoU A,_d+i = U - 1 - I,), 

• for n = 3, D ^3, X ^ AoU A,_^ U ^2^-4 = U (2z - 4 - I^) U A^_2 

• for the rest of cases X = U Az-d^\- 

4 The fractional part 

In this section will be shown the methodology to analyze - attractors of simple 
iterated functional systems. 
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4.1 The convex hull 

I will shortly present the results from [3], which allows to find analytically the 
convex hull of such simple fractals. 

The idea is to define some function which can be easily rewritten in the selfsimi- 
larity form - we get some functional equations, which can be solved analytically or 
approximated. 

The function we need in this case, is the width function, which gives for every 
direction the width in that direction. 

The other functions which can be threaded in such way are for example: 
/(x) = n {Fz + x)) (some Husdorff's measure now), 
fp{x) = T^yeiS^ - y)~^ (fo'^ example p = 2) 

/(^) = Ey6/z ^•''^■^(or = e^^'^dy) (this time '•' means scalar multiplication). 

In our cases, the equations for the width function can be solved analytically. 
They defines the convex hull of a set as the intersection of all halfplanes: 

conv(F, -Xo)= Pi {x : 3? (e-^°x) < h{a)} 

«e[0,27r) 



j>0 

In - 1 



Xo 



2z-l 

where xo is its center of symmetry, h{a) gives the position of bounding line 
in a direction. There is also shown how to construct analytically this con- 
vex hull from triangles, find the length of its boundary {2^/n + 2) or its area 

((^ - 1) Ei>o I sin(i(^)|y^"'y 

4.2 The construction of the boundary 

In this subsection, will be shown the methodology of using succeeding approxima- 
tions of Fz'. 

p _ E^=l,..,k z'^n + F, _ E.=i,....^ 

^^/c ^^^^ 
k>0 

the first form is true for any A; > - we can think about F^ as a discrete net of 
smaller copies of F^. 

The second: sum of growing family of such nets, shows that taking the limit 
/c ^ oo, we can 'forget' about tiles (copies of Fz) - approximations restores the 
whole set. 

We will introduce discrete versions of some topological properties, hke path. 
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boundary, connectiveness, show how to make the step to the next approximation 
and finally that in the limit, they really corresponds to their continuous equivalent. 

We will finally prove: 
Theorem 14. For periodic cases {z^ = Dz — n, D E 2(— -\/n, \/n)), we have: 

• if D ^ 0, we can split (intersecting in at most one point) 6Fz into six connected 
components - common parts with its 'neighboring tiles ', 

• Fz - connected, simply connected, 



• mt{F,) = F,. 

The purely imaginary case: z = j\fn (D = 0) can be easily solved ^ - even 
digits corresponds real part, odd: imaginary, F^ is just a rectangle. In this case 

+ Iz behave not like for the others: hexagonal lattice, but it's just a rectangular 
lattice - we will omit this case in this subsection. 

Look on ([9]) - changing the sign of D will only transpose F^ - in this subsection we 
can restrict to 

D <0. (16) 
For approximations we can work with X = 1^ zTL space: 

Definition 15. For any A C X define: 

Magnification of A is M{A) = Mz + n = Ui=o,..,n-i(^-2 + 0- 
k-th approximation of A is M^{A)/z^. 

We can think about it that we see only points from discrete lattice (X) and we 
use magnifying glass, which allows us to increase magnification \z\ times in one step 
(plus rotation). 

After such step, every point (tile) occurs to be n points(tiles). 



To have the hexagonal behavior, we would need to describe neighbors (Fig. [7|). 
I will just give positions and in the next lemma will be shown how to check them. 

Definition 16. 

Neighbors index set: R = Tjq (Z modulo 6) . 

r-th neighbor of a; G X (where r G i?) is Nr{x) := x + Nr, where: 

No ■= 0, Ni:=-z + D + 1, N2 := -z + D, N3 := -1, N^ := z - D ~ 1, N^ := z - D. 

The r-th meeting place is Cr{x) := xz + Cr, where 

Co := n + d, Ci:=n- 1, C2 := n - 1, C3 := -D - 1, C4 := 0, C5 := 0. 

The r-th meeting point is Pr{x) := x + Pj, where 

n-D r> . n tj . 1 id D-1 id. 1 



1 + ^,Pi := 1 + ^,P2 := ^,^3 := ^ - 1,P4 := ^ - 1,^5 

2—1' -'- 2—1 ' ^ 2—1' 2—1 ' * 2—1 ' "J 



2-1 ■ 
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Figure 7: The construction of magnification. 



N arc the positions of centers of neighboring tiles - we order them clockwise. 
After one step of magnification, the center tile is split into n tiles, Cr tells in which 
of them is the common point with the magnifications of r-th and r + 1-th neighbor 
(dots with numbers on the picture). 

We can easily find the real point of this meeting - do the magnification: 

zPr = Pr + the position of tile for which zPr is r-th place of meeting 

for example zPq = Pq + {n + D) + {z — D — 1). 

We will now check properties for making a step of magnification. 

Lemma 17. For x & X, r E R, we have that Cr{x) G M{x) is the neighbor of some 
elements from M{N^{x)) and M{Nr+i{x)). 

Proof: We will check it for x = - we can translate everything to get the thesis. 
It's easy to check the table below {z'^ — Dz — n): 
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Pr 
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1 + 

z—l 


n 
z-1 


^-1 


D 1 1 
z—l 


1 

z-1 



such that Nr^ — zN^ + a, A^^^ = zNr+i + h, 
for example for r = : Nr^(Cr) ^n + D + N4^z + (n-l)^ zNr + a □. 



Definition 18. 

The neighborhood of a; G X is N{x) := {Nr{x) : r € i?}. 
The ed^e of A C X is (UeA^(a^)) 

A sequence (xi)o<i<fc (xj e 5 C a;) will be called a length k chain in B if succeeding 
elements are neighbors. 

Chain is called closed, when additionally xq and x^-i are neighbors. 

Subchain of some closed, length k chain if is created from some of its succeeding 

elements (0 is after k — 1). 

A set A C X is called connected, if any two of its points can be connected by a chain 
in this set. 

A set A G X is called simply connected, if X\A is connected. 

The range of indexes: [a,b], where a,b E R, is called a step 1 arithmetic sequence 
in R from a to 6 (eg. [4, 2] = (4, 5, 0, 1,2)). 

Asei B <Z X lies on the left of sequence (a;j)o<i<fc, when for each 1 <i <k — 2: 
if Xi-i and Xj+i are correspondingly a and b neighbor of xi, then 

{xi + Nk : k e [a + l,b - 1]} G B. 

A set 5 C X lies on the left of closed chain, if this property is fulfilled for each of 
its subchains. 



We need lying on the left to make closed chain around some set oriented (coun- 
terclockwise) . 

The approximation of the boundary of can be started from length 6 closed chain 
- neighbors of in counterclockwise order - {0} lies on the left of this chain. 
Now we have to make a general construction of making the step to the next magni- 
fication. 

Lemma 19. Let B G X lies on the left of some (closed) chain (xj)o<i<A;- 

Then there exists (closed) chain {yi)o<i<i, such that yo G M{xo), G M{xk-i) 

and that P{B) lies on the left of it. 
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Proof: Firstly we will find for length tree chain: [Na, 0, A'^f,) and set 
{iV,+i,..,iVb_i}. 

Define i/q G M{Na) to be neighbor of Ca, yi-i G M{Ni,) to be neighbor of Cb- 
The rest of the chain lies in M(0) = n, such that the set lies on the left - it's chosen 
unambiguously - look on Fig. [7] : using some arrow we enter n in Ca, go along 
clockwise cycle inside, to exit turning left in Cb- 

For example (denote (iV„,0, A^f,) as Ofe): O4 ^ 4:o{3o)-^-%{0-i)''+^~^02. 

In the general situation, we are taking succeeding Xi and use the above construction 

for {xi-i — Xi, 0, Xj+i — Xi (remembering to remove duplicates on the ends). □ 

Now starting from Cq - closed, counterclockwise chain around {0}, call Cj - the 
z-th use of above lemma to Cq. 

It's the edge of M'({0}), which lies on the left of d. 

When we take approximations (divide this situation by z*), we see that we are 
getting closer to and its (oriented) boundary. 

C = {x e C : 3^^ Xi ^ X, XiE Ci/z'} 



Proof of Theorem I4 '■ We have only to check: 

• C is closed - just take diagonal sequence. 

• C is connected - the distance between succeeding points of chain Cj/z* goes 
to zero. 

• C\C has exactly two connected components - it's true for every approximation, 
using the previous point we have it in the limit. 

• The interior component with C is F^, C = 6Fz - M*({0})/z* are approxima- 
tions of Fz and its distance to Ci/z^ is going to 0. 

• In each step we can divide Ci/z^ into six sub chains - lying in approximations 
of succeeding neighbors {M\Na)/z'^) - we can split C into 6 subsets this way. 



• 'mi{Fz) = Fz - from the Bair's theorem: if int{z) = 0, then the sum of count- 
able number of them couldn't give the whole space. 

Because Vfc Fz = {M^{{0}) + Fz)/z^, we can cover Fz with its smaller copies 
having arbitrary small diameter - some of them will lie in this nonempty inte- 
rior. After rescaling back, we've got the thesis. □ 

Using the above construction of succeeding approximations of 5Fz, we can for 
example calculate its Hausdorff 's dimension. 

Namely using Lemma 19, we can change two edges: {Na,Q,Nb) into a sequence of 
them in the next magnification. 

If for each edge we distinguish between the direction of the next edge - remember 
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the first of tliis two edges as we can assign to it a sequence of such pairs in the 
next magnification (like in the example, but without the last O2). 
If we are interested in the number of edges only (instead of above grammar), we can 
write this iteration in the 36 x 36 (or 18 x 18 if we use symmetry) matrix form. It's 
normalized dominant eigenvector tells us the asymptotic probability distribution of 
edges, corresponding eigenvalue(A) tells that asymptotically, for the next approxi- 
mation we need A times more edges, but the length of them is y/n times smaller, so 
the boundary Hausdorff 's dimension is: 



H 



log(A) 
log(v^) 



Here are found values for n < 9: 



D/ 

n 





1 


2 


3 


4 


5 


2 


1 


1.210760533 


1.523627086 








3 


1 


1.162039854 


1.376841713 


1.657559542 






4 


1 


1.134761994 


1.303052340 


1.508664987 






5 


1 


1.116924317 


1.257583258 


1.422944863 


1.608726378 




6 


1 


1.104171451 


1.226294386 


1.366294523 


1.520716574 




7 


1 


1.094508825 


1.203216411 


1.325629733 


1.458928649 


1.598134771 


8 


1 


1.086882303 


1.185363533 


1.294784127 


1.412801706 


1.535582008 


9 


1 


1.080677473 


1.171064018 


1.270444950 


1.376841713 


1.487192945 



In our cases, 5Fz in fact - defines F^, and it's dimension is smaller. 
So we can use just the boundary 
faster. 



operating on it (drawing for example) is much 



5 Higher dimensions 

In this informal section there will be shortly introduce a suggestion of higher 

dimensions generalization, focusing on the dimension 3. 

Once more I will do it only for the periodic case. 

I haven't even checked the hypothesis numerically in this case. 



We can think about from the previous sections as 



z 



z 



where 1 = (1, 0)^, Z 

e this c 
We will focus on 



r 



cos{ip) — sm{(p) 
sin((y9) cos{ip) 



(17) 



We can use this definition (17) in higher dimension Euclidean space. 
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Choose some unit vector (l),rGM: |r|>l and some orthogonal matrix 

G M3x3_ 

Define 

Z := rO 

similarly as for complex numbers. 

We can diagonalize O - it is a rotation {ip radians) around some vector (u) and 
eventually the symmetry, but we can put it in the sign of r. 

We want 1, z, to generate whole space, so the angle between 1 and u (denote it 
-0) cannot be the multiplicity of 7r/2. 

So we can define orthonormal basis of our space: ei := 1, 62 is orthogonal to ei and 

u lies on the space generated by them, 63 is orthogonal to ei and 62- 

Define the another basis {u\v,w}: 

u' = cos{ip) ■ (cos('?/'), sin('?/'), 0) - rescalled u, 

V = sm{ip) ■ {sm{ip), — cosi^ip), 0) - orthogonal to u, 

w = 63. 

Now we can join v and w and think about our space as M x C: 

1 = (1, 1), we have multiplication and addition by coordinates, 

2^ = (r^ r'^e^'^^) (18) 

We can define the analog of the conjugancy: {61,62} plane symmetry 
(^{x,y) := {x,y)y 

Analogically to ([s]), we have the condition: 

N/ 

r = yn 

in N dimensional space. 

To assure the periodicity in dimension 3, we would need (|r|^ = n): 

= -Az^ -Bz-Cl (19) 

for some A,B,C G Z. 

conjugate this equation and multiply by z^: 

r^l = -Ar^z - Brh^ - Cz^ 
comparing both equations, we get (the sign of C is hidden in r): 

C = rye, A = Br^C, B = Ar^/C 
C = -r\ B = -Ar 
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We see that r have to be integer. We've found the analog of (13): 

(-r^-Ar,A,l), = (20) 

We now have to find some z = {R, Re^"^) fulfilhng this equation. 
The real part gives: 

= -r^ - ArR + AR^ + R^ 

ReR, so R = r. 

The complex part gives now: 

= -r^- Ar^e^^ + Ar^e^^^ + r'e^^^ = {r{e^^^ - 1) + Ae^'^{e^^ - 1)) 

-A/r = e-^''^(e2^'^ + e^''^ + 1) = e^''^ + 1 + e~^'^ = 1 + 2 cos(y?) 

(/? = arccos(^-^(l + A/r)^ A G {-3r + 1, r - 1}) (21) 

The sign of (f is only the matter of conjugation. 



Using (20) we can make sequence of reductions like in 3.1. 



By analogy to Lemma 1 1 , we have 



1^1 > -^|r(x)|<|x| (22) 

we will finally get into this ball, and finally get to a fixed point or a cycle. 
We can check that (0, 0, 0) is the only fixed point for A > 0. 
But in this case, for r > 0, we get a length two cycle: 
{2A - 1, -A. - 1) ^ {-A, -1, 0) ^ .... 
I've checked, that in the second case: 

m > 2, r = -m, n = m^, A G {0, 3m - 1}, (n, mA, A, 1)^ = (23) 

for a few first m everything is fine - we can generate the whole ball starting from 
(0, 0, 0) - we have proper numeral systems. 

We can generalize the methods for the fractional part in this cases too. 
The width function in introduced coordinates have two arguments, say 
a G [— 7r/2, 7r/2], /3 G [0,27r) which corresponds to (cos(a), sin(a)e-'^). 
Now using z'^Fz = + n(l, 1) + n{r,re^'^), we have equations for h, in which a is 
fixed - we can find solutions as a infinite sum for a fixed a, like in [3]. 
They don't correspond to one plane (cross section) this time - to find u coordinate, 
we need to use dh/da. 



We can use the method to construct the boundary of F^. 
This time every tile has 14 neighbors: 
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A^(0) = ±{1, A + z,A-l + z, Am + Az + z'^, Am -1 + Az + z"^, 

, Am - A + (A - l)z + z^, Am - A + 1 + {A - l)z + z'^} 
Now the simplest structure we can separately magnify, analogically to {Na,0, Nb), 
is with closed chain made of its neighbors, which splits N{0) into two connected 
subsets, which can be distinguished using chain orientation (one of them can be 
empty) . 
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